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Hydrostatic equilibrium of causally consistent and 
dynamically stable neutron star models 
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ABSTRACT 

We show that the mass-radius (M — R) relation corresponding to the stiffest equation 
of state (EOS) does not provide the necessary and sufficient condition of dynamical 
stability for the equilibrium configurations, since such configurations can not satisfy 
the 'compatibility criterion'. In this connection, we construct sequences composed of 
core-envelope models such that, like the stiffest EOS, each member of these sequences 
satisfy the extreme case of causality condition, v = c = 1, at the centre. We, thereafter, 
show that the M — R relation corresponding to the said core-envelope model sequences 
can provide the necessary and sufficient condition of dynamical stability only when 
the 'compatibility criterion' for these sequences is 'appropriately' satisfied. However, 
the fulfillment of 'compatibility criterion' can remain satisfied even when the M — R 
relation does not provide the necessary and sufficient condition of dynamical stability 
for the equilibrium configurations. In continuation to the results of previous study, 
these results explicitly show that the 'compatibility criterion' independently provides, 
in general, the necessary and sufficient condition of hydrostatic equilibrium for any 
regular sequence. Beside its fundamental feature, this study can also explain simul- 
taneously, both (the higher as well as lower) values of the glitch healing parameter 
observed for the Crab and the Vela-like pulsars respectively, on the basis of starquakc 
model of glitch generation. 

Key words: static spherical structures - dense matter - equation of state - stars: 
neutron - pulsars: individual: Crab: Vela. 



1 INTRODUCTION 

Characteristics of the super-dense objects like neutron stars 
(NSs) are based on the calculations of the EOS for the mat- 
ter at very high densities. However, the nuclear interactions 
beyond the density of ~ 10 14 gcm -3 are empirically not well 
known (Dolan f992) and the all known EOSs are only ex- 
trapolations of the empirical results far beyond this density 
range. In this regard, various EOSs based on theoretical ma- 
nipulations are available in the literature (Arnett & Bowers 
1977). Since the status of EOS of nuclear matter cannot be 
established empirically beyond a certain density range, one 
can apply physical constraints to obtain an upper bound of 
neutron star mass. Brecher and Caporaso (1976) assumed 
that the speed of sound in the nuclear matter equals the 
speed of light and obtained a value of 4.8M© as an up- 
per limit of the neutron star masses. However, the matter 
described by this stiffest EOS, (dP/dE) = 1 (geometrized 
units) or P — (E — E s ) [where P is the pressure, E is the 
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energy-density and E s — 2 x 10 1 gem -3 represents the sur- 
face density] has a super-dense self-bound state at the sur- 
face of the configuration where the pressure is vanishingly 
small, which represent the 'abnormal state of matter' (Lee 
1975; Haensel & Zdunik 1989). This 'abnormality' may be 
specified as - (i) the pressure vanishes at the average nuclear 
densities, and (ii) the speed of sound is equal to that of light 
even when the pressure is vanishing small. This 'abnormal- 
ity' can be removed if we ensure continuity of pressure, den- 
sity and both of the metric parameters at the boundary of 
the structure (Negi & Durgapal 2000). 

Earlier, Rhoades and Ruffini (1974), without going into 
the details of the nuclear interactions, assumed that beyond 
a certain density 4.6 x 10 14 gcm -3 (the range of densities 
where no extrapolated EOS is known) the EOS in the core 
is given by the criterion that the speed of sound attains the 
speed of light, that is, (dP/dE) = 1, and matched the core 
to an envelope with the BPS (Baym, Pethick & Sutherland 
1971) EOS and obtained an upper limit for the neutron star 
mass as 3.2Mq . Hartle (1978) emphasized that the max- 
imum masses of neutron stars obtained in this manner in- 
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volve a scale factor, k — [_E m /10 14 gcm _3 ] _1 ' /2 , such that, 
the matching density, E m , plays a sensitive role to obtain 
an upper bound on neutron star masses. Usually, k is taken 
to be equal to or greater than one in all the conventional 
models. For the densities less than E m , the matter compos- 
ing the object is assumed to be known and unique. That is, 
the EOS of the envelope of these stars are chosen so that the 
'abnormalities' in the sense mentioned above are removed. 
Friedman and Ipser (1987) calculated the masses of the neu- 
tron stars for different values of the matching densities using 
EOS given by BPS and NV (Negele & Vautherin 1973), re- 
spectively, and concluded that the EOS chosen for the enve- 
lope does not make any significant difference in the results, 
because in each case the mass in the envelope turns out to 
be insignificant as compared to that of the core containing 
the most stiff material. 

In order to implement this 'insignificant' mass of the 
envelope component, Crawford & Demiahski (2003) have re- 
cently constructed the NS models for seven representative 
EOSs of dense nuclear matter by covering a range of NS 
masses. They have computed the 'fractional moment of in- 
ertia' of the core component which is defined as the glitch 
healing parameter, Q, in the starquake mechanism of glitch 
generation as 

Q=^I±, (1) 

J total 

where /total represents the moment of inertia of the entire 
configuration. Their study shows that the much larger val- 
ues of Q(> 0.7) for the Crab pulsar are fulfilled by all but 
the six EOSs considered in the study corresponding to a 're- 
alistic' neutron star mass range 1.4 ± O.2M . On the other 
hand, for the much lower values of Q(< 0.2) correspond- 
ing to the case of the Vela pulsar, their models predict a NS 
mass < 0.5M & which is too low as compared to the 'realistic' 
NS mass range. Thus, their study concludes the 'starquake' 
as a feasible mechanism for glitch generation in the Crab- 
like pulsars (corresponding to a larger value of Q) and the 
'vortex unpinning', the another mechanism of glitch gener- 
ation, suitable for the Vela-like pulsars (corresponding to a 
lower values of Q). However, it seems quite surprising that 
if the internal structure of NSs is supposed to be described 
by the same two-component conventional model, then why 
different kinds of glitch mechanisms are required to explain 
a glitch! 

This is the purpose of this study to provide an insight 
about a fundamental 'theorem' concerning the hydrostatic 
equilibrium of NS models (which was laking in the earlier 
studies) and to show that as soon as this theorem is imple- 
mented to an appropriate conventional NS sequence, the var- 
ious shortcomings of the conventional NS models (as men- 
tioned above) can be resolved. 



2 REMOVAL OF ABNORMALITIES FROM 
THE STIFFEST EOS AND BOUNDARY 
CONDITIONS FOR NEUTRON STAR 
MODELS COMPATIBLE WITH 
HYDROSTATIC EQUILIBRIUM 

In order to construct such an appropriate sequence of NS 
models, consistent with causality and dynamical stability, 



we offer here an entirely different approach to the whole 
problem which will not only remove the 'abnormalities' of 
the stiffest EOS (as discussed earlier in detail under section 
1) but can also assure the necessary and sufficient condi- 
tion of hydrostatic equilibrium for the resulting configura- 
tion. Since, as we would show later in section 3 that the 
M — R relations corresponding to the configurations (i) gov- 
erned by the pure stiffest EOS, and (ii) those resulting from 
the removal of 'abnormalities' from the stiffest EOS (core- 
envelope models), do not provide the necessary and suffi- 
cient condition of dynamical stability unless the 'compati- 
bility criterion' (Negi & Durgapal 2001) which states that: 
"for each and every assigned value of a[= (Po/Eo) = the 
ratio of central pressure to central energy-density] , the com- 
pactness ratio u(= M/R) of the entire configuration should 
not exceed the compactness ratio, Uh, of the corresponding 
homogeneous density sphere (that is, u < uh)" is 'appropri- 
ately' satisfied by such configurations. 

The present approach is based on the 'theorem' applica- 
ble to a wide range of conventional NS sequences (including 
the core-envelope models) provided that every member of 
such sequences satisfies the condition (dP/dE)o = 1 (here 
and elsewhere in the paper, the subscript '0' represents the 
value of the corresponding quantity at the centre of the con- 
figuration). This theorem asserts that in order to assure the 
necessary and sufficient condition of dynamical stability for 
the mass, the maximum stable value of u (corresponding to 
the case of first maxima among masses in the mass-radius 
(M — R) relation) and the corresponding central value of 
'local' adiabatic index (ri) (= [{P + E )/P ](dP/dE) ) 
of such sequences must satisfy the inequalities, it max < 
iWx.abs = 0.3406 and (ri)o < (ri)o,max,abs = 2.5946 re- 
spectively (Negi 2007). In addition to the result of previous 
study regarding the 'compatibility criterion', we would show 
in the present study that "the M — R relation (or the mass- 
central density relation) corresponding to an equilibrium se- 
quence provides the necessary and sufficient condition of dy- 
namical stability only when the 'compatibility criterion' for 
the equilibrium sequence is also satisfied. But the fulfillment 
of the 'compatibility criterion' does not depend upon the ful- 
fillment of the necessary and sufficient condition of dynam- 
ical stability provided by the M — R relation" . The proof of 
this statement would explicitly show that the fulfillment of 
'compatibility criterion' independently provides, in general, 
the necessary and sufficient condition of hydrostatic equilib- 
rium for any sequence composed of regular configurations, 
since the first sentence of this statement has already been 
proved in the previous study for the equilibrium sequences 
of the type mentioned above (Negi 2007) . 

In order to provide a proof of the last statement, the 
present paper deals with the construction of a four step con- 
secutive study in the following manner 

Step (1): We would construct the NS models corre- 
sponding to the pure stiffest EOS, dP/dE = 1, and show 
that the M — R relation corresponding to this sequence does 
not provide the necessary and sufficient condition of dynam- 
ical stability, since the 'compatibility criterion' can not be 
fulfilled by the equilibrium configurations corresponding to 
this 'abnormal' EOS. 

Step (2): We would construct the NS sequences com- 
posed of core-envelope models by considering the stiffest 
EOS, dP/dE = 1, in the core and the EOS of classical 
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polytrope, dlnP/dlnp = Ti (where p is the density of the 
rest-mass and Ti is the constant adiabatic index; see, e.g. 
Tooper 1965), in the envelope (so that each member of these 
sequences satisfy the extreme case of causality condition, 
v = c = 1, at the centre) for an 'arbitrarily' assigned value of 
pressure to density ratio, (P b /E b ) = (P b /E b )i (say), at the 
core-envelope boundary. Though this procedure removes the 
'abnormalities' from the stiffest EOS, but as we would show 
later in sec. 3 that the M — R relation corresponding to the 
NS models with an envelope Fi = 5/3 and 2, do not provide 
the necessary and sufficient condition of dynamical stabil- 
ity, since the 'compatibility criterion' remains unsatisfied by 
various stable (as well as unstable) models corresponding to 
the said sequences. 

Step (3): In this sub-section we would re-construct the 
NS sequence composed of core-envelope models as described 
in step (2) for the same boundary value of (P b /E b )\. But in- 
stead of the constant Ti = 5/3 and 2, the constant Ti = 4/3 
would be used for the envelope. We would show that al- 
though the 'compatibility criterion' is satisfied by all stable 
(as well as unstable) models comprising the sequence but 
still, the M — R relation does not provide the necessary and 
sufficient condition of dynamical stability, since the 'com- 
patibility criterion' is not 'appropriately' satisfied by the 
equilibrium configurations. 

Step (4): We would re-construct the NS sequences de- 
scribed in steps (2) and (3) respectively for an 'appropri- 
ate' boundary value of (P b /E b ) — (P b /E b )2 (say). Since 
for this particular value of (P b /E b ), The 'compatibility cri- 
terion' turns out to be 'appropriately' satisfied by all the 
equilibrium sequences corresponding to an envelope with 
Ti = 4/3,5/3 and 2 respectively, it would follow that the 
M — R relation does provide the necessary and sufficient con- 
dition of dynamical stability for the said sequences [the value 
(P b /E b )'2 corresponds to the minimum value of (P b /E b ) for 
which the M — R relation fulfills the necessary and sufficient 
condition of dynamical stability for the equilibrium sequence 
corresponding to NS models with an envelope I\ = 2. Since 
the value (P b /E b )2 turns out to be greater than (P b /E b )i, it 
follows, therefore, that this particular value (P b /E b )2 would 
guarantee the automatic fulfillment of the 'compatibility cri- 
terion' for all the sequences corresponding to NS models 
with an envelope Ti = 4/3, 5/3 and 2 respectively, together 
with fulfilling the necessary and sufficient condition of dy- 
namical stability provided by the M — R relation. This con- 
dition is termed as the 'appropriate' fulfillment of 'compat- 
ibility criterion'.]. 

The methodology regarding the construction of various 
models mentioned in steps (1) - (4) above and the important 
outcome emerging thereafter are discussed in the following 
section. 

The assignment of different values of Fi between the 
range (4/3) and 2 in the envelope of NS models discussed un- 
der step (2) to step (4) above follows from the fact that (4/3) 
represents the EOS of ultrarelativistic degenerate electrons 
and non-relativistic nuclei (Chandrasekhar 1935) or of rela- 
tivistic degenerate neutron gas (Zeldovich & Novikov 1978). 
(5/3) corresponds to the well known EOS of non-relativistic 
degenerate neutrons gas (Oppenheimer & Volkoff 1939) and 
Fi — 2 represents the case of extreme relativistic baryons in- 
teracting through vector meson field (Zeldovich 1962) [The 
value of Fi > 2 is also possible for some EOS, e.g., Mal- 



one, Johnson & Bethe (1975); Clark, Heintzmann & Grow- 
ing (1971), however, the results obtained in this paper do 
not get affected by choosing Ti > 2]. The various values of 
Ti chosen in this range (4/3 < Fi < 2) can cover almost all 
of the nuclear EOS discussed in the literature (and, there- 
fore, it can cover almost full range of nuclear density which 
might be applicable for the envelope region as well), and in 
our opinion, it could be more appropriate to choose an av- 
erage (constant) value of Ti for a conventional NS model, 
instead of going into the details of the density range be- 
low E b specified by different EOS [for example, BPS, NV, 
or FPS (Lorenz, Ravenhall & Pethick 1993)] which are fre- 
quently used by various authors in the conventional mod- 
els of neutron stars in-spite of their uncertainty (Friedman 
& Ipser 1987; Dolan 1992) as mentioned earlier. The EOS 
of classical polytrope for the said Ti values, considered in 
the present study, will not only simplify the procedure but 
also provide necessary insight regarding the suitability of 
the EOS for the envelope region, since this study yields the 
important finding that the stable sequences of neutron star 
models terminate at the same value of maximum mass in- 
dependent of the EOS of the envelope. 



3 METHODOLOGY AND DISCUSSION 

The metric for spherically symmetric and static configura- 
tions can be written in the following form (remembering that 
we are using geometrized units, i.e. G = c = 1; where G and 
c represent respectively, the universal constant of gravitation 
and the speed of light in vacuum) 



e at 



(2) 



where v and A are functions of r alone. The Oppenheimer- 
Volkoff (O-V) equations (Oppenheimer & Volkoff 1939), 
resulting from Einstein's field equations, for systems with 
isotropic pressure P and energy-density E can be written as 



P 1 = - (P + E) [4irPr 3 + m]/r(r- 
v'/2 = -P'/(P + E) 
m'(r) = AnEr 2 ; 



2m) 



(3) 
(4) 
(5) 



where m(r) = J" 7 " 4ivEr 2 dr is the mass, contained within the 
radius r, and the prime denotes radial derivative. 

In order to solve equations (3) - (5) for different models 
considered in the present study, we consider the stiffest EOS 



P = E-E S 



(6) 



for the entire configuration (or for the core region), and the 
EOS of classical polytrope 

dlnP 



dlnp 



= Ti 



(7) 



for the envelope region respectively. For core-envelope mod- 
els, the fractional moment of inertia (given by equation (1)) 
may be calculated by using an approximate but very precise 
empirical formula which is based on 30 theoretical EOSs of 
dense nuclear matter. For NS models, the formula yields in 
the following form (Bejger & Ffaensel 2002) 



I~ -(l + 5x)MR 2 , x>0.l 



(8) 
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Figure 1. The mass-radius diagram of the models correspond- 
ing to the stiffcst EOS, P = (E — E s ), (as discussed in the text 
(step 1, sections 2 and 3)) for an assigned value of the density 
E 3 = 2.7 X 10 14 g cm~ 3 at the surface. The models do not ful- 
fill the 'compatibility criterion' as shown in Table 1. Also, the 
mass-radius diagram does not provide the necessary and suffi- 
cient condition of dynamical stability because the inequalities, 
(Fi)o < 2.5946 and ti max < 0.3406, are not fulfilled simultane- 
ously at the maximum value of mass as shown in Table f . 

where x is the compactness parameter measured in units 
of [M (km)/km], i.e. 



M/R u 
M Q /km ~~ 1.477' 



(9) 



In accordance with the four step consecutive study men- 
tioned in the last section, we construct different models in 
the following manner 

Models mentioned in step (1): The coupled equa- 
tions (3 - 5) are solved for the EOS given by equation (6) 
until pressure vanishes at the surface of the configuration for 
a fiduciary choice of E s = 2.7 x 10 14 g cm -3 , the nuclear sat- 
uration density. At the surface, r = R, we obtain P = 0, E = 
E s ,m(r = R) = M,e" = e~ x = (1 - 2M/R) = (1 - 2u). 
The results of the calculations are shown in Table 1 and the 
M — R diagram is presented in Fig.l. It follows from Table 
1 that along the stable branch of the sequence, the max- 
imum value of mass corresponds to the maximum 'stable' 
value of u max — 0.3539 and the corresponding 'local' value 
of (Ti)o ~ 2.4990. Although, this value of (Ti) turns out 
to be consistent with that of the absolute upper bound on 
(ri)o((Ti)o, m ax,abs — 2.5946), the maximum 'stable' value of 
u max ~ 0.3539 is found to be inconsistent with that of the 
absolute Upper bound OU Wmax(Wmax,abs =< 0.3406). Thus the 
configuration turns out to be inconsistent with the corollary 
1 of Theorem 2 (Negi 2007). It follows, therefore, that the 
M — R relation corresponding to the stiffest EOS does not 



Figure 2. The mass-radius diagram of the models as discussed 
in the text (steps (2) and (3), sections 2 and 3) for an assigned 
value of matching density E = E b = 2.7 X 10 14 g cm~ 3 at 
the core-envelope boundary. The labels a, b and c represent 
the models for an envelope with IT = (4/3), (5/3), and 2 re- 
spectively. The value of the ratio of pressure to energy-density, 
(P b /E b ) = (P b /E b ) 1 = 1.0645X10" 2 , at the core envelope bound- 
ary is assigned in such a manner that the 'compatibility criterion' 
has become satisfied by all the models corresponding to an en- 
velope with IT = 4/3, whereas it remains unsatisfied by various 
models corresponding to an envelope with Ti =5/3 and 2 respec- 
tively as shown in Table 3. Since, the 'compatibility criterion' has 
not been satisfied 'appropriately' by the models corresponding to 
an envelope with IT = 4/3 (and has not been satisfied by various 
models corresponding to an envelope with Ti = 5/3 and 2), the 
M — R relation does not provide the necessary and sufficient con- 
dition of dynamical stability for any sequence as shown in Table 
2. 



provide the necessary and sufficient condition of dynamical 
stability for the equilibrium configurations. Since, the total 
mass 'M' which appears here does not fulfill the definition 
of the 'actual mass' which should be present in the exterior 
Schwarzschild solution (Negi 2004; 2006), the equilibrium 
sequence corresponding to the pure stiffest EOS, therefore, 
does not even fulfill the necessary condition of hydrostatic 
equilibrium (Negi 2007), as a result the 'compatibility crite- 
rion' can not be made satisfied by such configurations. This 
is also evident from the comparison of column 2 and 6 of 
Table 1 that for each assigned value of (Po/Eo), u s tg > Uh- 
Models mentioned in step (2): The coupled equa- 
tions (3 - 5) are solved by considering equation (6) in the core 
(0 < r < b) and equation (7) in the envelope (b < r < R) 
for the constant Ti = 5/3 and 2 respectively, until the 
boundary conditions: P — E = 0, m(r — R) = M, e v — 
e~ A = (1 — 2M/R) = (1 — 2u) are reached at the surface, 
r = R, of the configuration. Equation (8) is also used to- 
gether with equations (3 - 5) to calculate the fractional mo- 
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Figure 3. Fractional moment of inertia Q(= ^core/Aotal) vs - to- 
tal mass M for the models presented in Tables 2 and 3 and shown 
in Figure 2. The labels a, b and c denote the models for an enve- 
lope with Ti = (4/3), (5/3), and 2 respectively. 



ment of inertia given by equation (1) and the moment of 
inertia of the entire core-envelope model considered in the 
present section. The calculations are performed for a fidu- 
ciary choice of the boundary density, Eb — 2.7 x 10 14 g 
cm -3 , and for an 'arbitrarily' assigned value of the ratio of 
pressure to energy-density at the core-envelope boundary, 
(P b /E b ) = (P b /E b )i = 1.0645 x 1(T 2 . The results of the cal- 
culations are summarized in Tables 2-3 (indicated with a 
superscript 'b' for the Ti = 5/3 and 'c' for the Fi = 2 enve- 
lope models) and the M — R diagram is shown in Fig. 2 (label 
'b' for the Fi = 5/3 and V for the Fi = 2 envelope models). 
Though, the construction of such core-envelope models ap- 
propriately fulfill the definition of the 'actual mass' appear- 
ing in the exterior Schwarzschild solution (Negi 2007), the 
choice of (Pb/E b ) = (Pb/E b )i, however, yields the sequences 
(corresponding to NS models with an envelope Ti = 5/3 and 
2 respectively) which fulfill only the necessary (but not suf- 
ficient) condition of hydrostatic equilibrium, since the 'com- 
patibility criterion' remains unsatisfied by such sequences 
as shown in Table 3. From the Table 2, we find that at 
the maximum mass values along the stable branch of the 
sequences, the sequence composed of NS models with an en- 
velope Fi = 5/3 yields Umax — 0.3493 and the corresponding 
value of (ri)o — 2.4984, whereas the sequence composed of 
NS models with an envelope Ti = 2 yields u max ~ 0.3509 
and the corresponding value of (Fi)o ~ 2.4984 respectively. 
Both pairs of these values are, however, found to be in- 
consistent with that of the pair of absolute upper bounds 
Umax,abs(= 0.3406) and (ri) 0jmaXiabs (^ 2.5946) respectively. 
It follows, therefore, that the M — R relation corresponding 
to the sequences (composed of NS models with an envelope 



Ti = 5/3 and 2 respectively) does not provide the necessary 
and sufficient condition of dynamical stability as shown in 
Fig.2. 

Tables 2-3 and Fig. 3 show that the range of frac- 
tional moment of inertia, 0.652 < Q < 0.948, correspond- 
ing to the 'stable' models (with an envelope Ti = 5/3 
and 2 respectively) possesses the NS masses in the range 
1.35M n < M < 4.1M ffl . This feature is consistent with 



those of the conventional models discussed in the literature 
and can explain only the higher values of the glitch healing 
parameter observed for the Crab-like pulsars. However, for 
the minimum weighted mean value of Q ~ 0.7 correspond- 
ing to the Crab pulsar, we obtain from Fig. 3, the minimum 
masses M b ~ 1.6M and M c ~ 1.44M© for the Crab pulsar. 

Models mentioned in step (3): In order to con- 
struct core-envelope models corresponding to the constant 
Ti = 4/3 for the envelope, we use equation (7) for the 
envelope (6 < r < R) and equation (6) for the core 
(0 < r < b) and solve equations (3 - 5) together with 
equation (8) for the same boundary conditions mentioned 
in the last sub-section (step (2)) for the surface and the 
core-envelope boundary, viz. E b = 2.7 x 10 14 g cm~ 3 , and 
(Pb/E b ) = (P b /E b )i = 1.0645 x 10~ 2 respectively. The vari- 
ous parameters obtained for this model are indicated with a 
superscript 'a' in Tables 2 and 3 respectively, and the M — R 
diagram (marked with label 'a') is presented in Fig.2. Table 3 
shows that although the 'compatibility criterion' is satisfied 
by all members of the sequence corresponding to NS models 
with an envelope Fi = 4/3 for the choice of the bound- 
ary condition (P b /E b ) = (Pb/E b )\, this choice, however, 
does not appear to be 'appropriate' because at the maxi- 
mum value of mass along the stable branch of the sequence 
in the mass-radius relation, the model yields the maximum 
value of u(u max ) — 0.3444 and the corresponding value of 
(ri)o ~ 2.4984 as shown in Table 2. These values, however, 
show inconsistency with that of the pair of absolute upper 
bounds, u max ,abs(= 0.3406) and (ri)o, max , abs (= 2.5946), re- 
spectively. It follows from this result that the necessary and 
sufficient condition of dynamical stability (provided by the 
M — R relation) may remain unsatisfied even when all mem- 
bers of the sequence do satisfy the 'compatibility criterion'. 

Tables 2-3 and Fig.3 also indicate that the range 
of fractional moment of inertia, 0.558 < Q < 0.922, for 
the 'stable' sequence corresponding to the Y\ = 4/3 enve- 
lope model provides the masses of NS models in the range, 
1.36M© < M < 4.1M Q . These higher values of the glitch 
healing parameter are consistent only with those the Crab- 
like pulsars as discussed in the last sub-section (step 2). 
However, for the minimum value of Q ~ 0.7, Fig.3 yields 
the minimum mass M c ~ 2.O3M0 for the Crab pulsar. 

Models mentioned in step (4): We re-construct the 
NS sequences (composed of core-envelope models), those 
described in the last two sub-sections (steps (2) and (3) 
respectively), but for an 'appropriate' value of (Pb/Eb) = 
(Pb/E b ) 2 = 4.694 x 10~ 2 . The results of the calculations 
are summarized in Tables 4-5 and the M — R diagram is 
shown in Fig. 4. The superscripts 'a', 'b' and 'c' which ap- 
pear among various parameters in these tables, represent the 
models with an envelope I\ = 4/3,5/3 and 2 respectively. 
The M — R diagram is also labeled in the similar fashion. 
The value (P b /E b )2, in fact, represents the minimum value 
of (Pb/Eb) for which all the sequences, corresponding to NS 
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Figure 4. The mass-radius diagram of the models as discussed 
in the text (step 4, sections 2 and 3) for an assigned value of 
matching density E = Eb = 2.7 X 10 14 g cm -3 at the core- 
envelope boundary. The labels a, b and c represent the mod- 
els for an envelope with Ti = (4/3), (5/3), and 2 respectively. 
The value of the ratio of pressure to energy-density, (P b /Eb) = 
(Pb/E b ) 2 = 4.694 x 10" 2 , at the core envelope boundary is ob- 
tained in such a manner that the 'compatibility criterion' is 'ap- 
propriately' satisfied by all the models corresponding to an en- 
velope with Ti = 4/3, 5/3 and 2 respectively, as shown in Table 
5. That is, the necessary and sufficient condition of dynamical 
stability provided by the M — R relation is also satisfied together 
with satisfying the 'compatibility criterion' as shown in Table 4. 



models with an envelope T\ = 4/3,5/3 and 2 respectively, 
fulfill the necessary and sufficient condition of dynamical 
stability provided by the M — R relation. Since the value 
(P b /Eb)2 > (Pb/Eb)i, it follows that the 'compatibility cri- 
terion' would automatically be satisfied by all the sequences 
corresponding to NS models with an envelope 1\ = 4/3, 5/3 
and 2 respectively. This is evident from the Tables 4 and 5 
respectively. A comparison of column 2 with those of col- 
umn 3, 4 and 5 of Table 5 show that the 'compatibility cri- 
terion' is 'appropriately' satisfied by all the sequences. Table 
4 shows that at the maximum values of mass along the sta- 
ble branch of the sequences, the models corresponding to 
an envelope with Fi = 4/3, 5/3 and 2 yield the maximum 
values of w(w mas ) — 0.3044, 0.3324 and 0.3404 respectively 
for the same value of (Ti)o — 2.4955 . Evidently, these val- 
ues are consistent with those of the absolute upper bounds 
0.3406) and (Ti) 

0,max,abs 

(2* 2.5946) respectively. 
Tables 4-5 and Fig. 5 show that the ranges of fractional 
moment of inertia, 0.069 < Q < 0.779, obtained for the sta- 
ble sequences corresponding to the Ti = 5/3 and 2 envelope 
models have masses in the range, 1.48M0 < M < 4.2Mq. 
Obviously, these ranges are capable of fulfilling the higher 
as well as lower values of the glitch healing parameter cor- 



Figure 5. Fractional moment of inertia Q(= ^core/Aotal) vs - to- 
tal mass M for the models presented in Tables 4 and 5 and shown 
in Figure 4. The labels a, b and c denote the models for an enve- 
lope with Ti = (4/3), (5/3), and 2 respectively. 



responding to the case of the Crab and the Vela-like pulsars 
respectively. For the maximum weighted mean value of Q ~ 
0.2 corresponding to the Vela pulsar, Fig. 5 yields the maxi- 
mum masses M b ~ 2.06M Q and M c ~ 1.85M for the Vela 
pulsar, whereas for the minimum weighted mean value of 
Q ~ 0.7, the models yield the minimum masses Mb — A.IMq 
and M c ~ 3.98M0 for the Crab pulsar respectively. It also 
follows from the Tables 4-5 and Fig. 5 that for the stable 
range of fractional moment of inertia, 0.016 < Q < 0.655, 
corresponding to the Y\ = 4/3 envelope models, the NSs can 
have masses in the range, 1.79M Q < M < 4.2M Q . For the 
maximum weighted mean value of Q ~ 0.2, Fig. 5 yields the 
maximum mass M c ~ 2.67M© for the Vela pulsar, however 
the minimum weighted mean value, Q ~ 0.7, corresponds 
to the minimum mass, M a w 3.98M0, for the Crab pulsar 
which lies in the unstable branch of Fig. 5. 

The important conclusions emerge from the four step 
study and discussion presented above may be summarized 
in the following manner: 

(1) Steps (1), (2) and (4) show that the M — R rela- 
tion provides the necessary and sufficient condition of dy- 
namical stability only when the 'compatibility criterion' for 
the equilibrium configurations is satisfied. Without satisfy- 
ing the 'compatibility criterion', the M — R relation can not 
provide the necessary and sufficient condition of dynamical 
stability. 

(2) The inclusion of step (3), however, modifies the 
above conclusion in the following statement: The M — R 
relation provides the necessary and sufficient condition of 
dynamical stability only when the 'compatibility criterion' 
is 'appropriately' satisfied. Thus, the necessary and suffi- 
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1.5 

(M/M e ) 



Figure 6. Fractional moment of inertia Q(= /core //total ) vs - 
total mass M for the models as discussed in the text (step 
4, section 2 and 3) for a calculated value of matching density 
E = Ei, = 7.0794 X 10 14 g cm -3 at the core-envelope boundary. 
The labels a, b and c represent the models for an envelope with 
Ti = (4/3), (5/3), and 2 respectively. The value of the ratio of 
pressure to energy-density, (P b /E b ) = (P b /E b ) 2 = 4.694 X 10~ 2 , 
at the core envelope boundary is obtained in such a manner that 
the 'compatibility criterion' is 'appropriately' satisfied for all the 
models corresponding to an envelope with Ti = 4/3,5/3 and 2 
respectively (in the same manner as shown in Table 5 for the 
matching density E = E b = 2.7 X 10 14 g cm -3 at the core- 
envelope boundary). That is, the necessary and sufficient condi- 
tion of dynamical stability provided by the M — R relation is also 
satisfied together with satisfying the 'compatibility criterion' (in 
the same manner as shown in Table 4 for the matching density 



E = E b = 2.7 X 10 g cm 3 at the core-envelope boundary). 



cient condition provided by the M — R relation is totally 
dependent on the 'appropriate' fulfillment of 'compatibil- 
ity criterion'. On the other hand, the fulfillment of 'com- 
patibility criterion' is quite independent of the fulfillment 
of necessary and sufficient condition of dynamical stability 
provided by the M — R relation, because the 'compatibil- 
ity criterion' can be made satisfied even without fulfilling 
the necessary and sufficient condition of dynamical stabil- 
ity provided by the M — R relation. Since, the fulfillment 
of the 'compatibility criterion' alone is a measure of the ful- 
fillment of necessary and sufficient condition of hydrostatic 
equilibrium for any static and spherical configuration (Negi 
2007), it follows from the above discussion that the fulfill- 
ment of the 'compatibility criterion' alone, independently 
provides, in general, the necessary and sufficient condition 
of hydrostatic equilibrium for any sequence composed of reg- 
ular, static and spherical configurations. 



We have investigated that the mass-radius relation corre- 
sponding to the stiffest equation of state does not provide 
the necessary and sufficient condition of dynamical stabil- 
ity because at the maximum value of mass along the mass- 
radius relation, the pair of the maximum 'stable' value of 
compactness, tt m ax — 0.3539, and the corresponding central 
value of the 'local' adiabatic index, (ri)o ^ 2.4990, turns 
out to be inconsistent with that of the pair of absolute val- 
ues, u m ax,abs = 0.3406, and (r x ) ,max,abs = 2.5946, compat- 
ible with the structure of general relativity, causality, and 
dynamical stability. The reason behind this inconsistency 
lies in the fact that the 'compatibility criterion' (Negi and 
Durgapal 2001) can not be fulfilled by any of the sequence, 
composed of regular configurations corresponding to a sin- 
gle EOS with finite (non-zero) values of surface and central 
density (Negi 2004; 2006). 

We have further constructed the sequences composed 
of core-envelope models such that, like the stiffest EOS, 
each member of these sequences satisfy the extreme case 
of causality condition, v = c = 1, at the centre by consider- 
ing the stiffest equation of state in the core and a poly- 
tropic equation with the constant adiabatic index Ti = 
[AlnP/Alnp] in the envelope and investigated the result- 
ing configurations on the basis of the 'compatibility cri- 
terion' and the mass-radius relation for different values of 
Ti = 4/3,5/3 and 2. Together with the finding correspond- 
ing to the case of stiffest EOS mentioned above, the inves- 
tigation of the said sequences of the core-envelope models 
explicitly show that the mass-radius relation corresponding 
to the said sequences can provide the necessary and sufficient 
condition of dynamical stability only when the 'compatibil- 
ity criterion' for the sequences in 'appropriately' satisfied. 
However, the fulfillment of 'compatibility criterion' can re- 
main satisfied even when the mass-radius relation does not 
provide the necessary and sufficient condition of dynamical 
stability for the equilibrium configurations. 

In continuation to the results of previous study (Negi 
2007), these results explicitly show that the 'compatibility 
criterion' independently provides, in general, the necessary 
and sufficient condition of hydrostatic equilibrium for any 
regular sequence. 

The causal sequences of NS models constructed in the 
present study (which fulfill the necessary and sufficient con- 
dition of hydrostatic equilibrium and dynamical stability 
simultaneously) terminate at the same value of maximum 
mass, M max « 4.2M (for E b = 2.7 x 10 14 g cm -3 ), inde- 
pendent of the EOS of the envelope. However, the maximum 
compactness, u max — 0.3404, yields for the sequence corre- 
sponding to the Y\ = 2 envelope models. While for the the 
same value of transition density E b = 2.7 x 10 14 g cm~ 3 at 
the core-envelope boundary, this upper bound on NS mass 
found to be fully consistent with those of the models formu- 
lated with the advance nuclear theory (Kalogera & Baym 
1996), the Y\ — 2 envelope model indicates the appropriate- 
ness of the (average) value of Fi for the entire envelope. 

Beside its fundamental feature, this study also under- 
lines the importance of the applicability of 'compatibility 
criterion' to the conventional models of NSs. Since, we find 
that when the 'compatibility criterion' is not satisfied (for 
the case of the models corresponding to an envelope with 
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T\ — 5/3 and 2, if the ratio of pressure to energy-density 
at the core-envelope boundary, Pb/Eb, is 'arbitrarily' as- 
signed to be about 1.0645 x 10~ 2 ) or not 'appropriately' 
satisfied (for the case of the models corresponding to an en- 
velope with Ti = 4/3, if P b /E b ~ 1.0645 x 10~ 2 ) by the 
sequences, the corresponding range of the glitch healing pa- 
rameter turns out to be 0.558 < Q < 0.948. This feature 
is consistent with the other conventional NS models dis- 
cussed in the literature and can explain only the higher val- 
ues of Q on the basis of starquake model of glitch generation 
for the Crab-like pulsars. For the minimum weighted mean 
value of Q ~ 0.7, our models (corresponding to an envelope 
with Ti = 2,5/3 and 4/3 respectively) yield the minimum 
masses M c ~ 1.44M©, M b ~ 1.6M© and M a ~ 2.O3M for 
the Crab pulsar. Among these values, the first two values 
are comparable with those of the minimum values 1.35Mq 
and 1.65Mq obtained by Crawford and Demiahski (2003) 
by using the GWM (Glendenning, Weber, & Moszkowski 
1992) and HKP (Haensel, Kutschera, & Proszynski 1981) 
EOSs. Whereas, corresponding to Q ~ 0.7, the other five 
EOSs of the dense nuclear matter considered by Crawford 
and Demiahski (2003) yield the minimum mass M < IMq 
for the Crab pulsar (see, e.g. Crawford and Demiahski 2003; 
and references therein) . However for the maximum weighted 
mean value of Q ~ 0.2 corresponding to the Vela pulsar, our 
models also yield the unrealistically small mass values for the 
Vela pulsar together with the other models discussed in the 
literature (Crawford and Demiahski 2003). 

On the other hand, as soon as the 'compatibility cri- 
terion' is 'appropriately' satisfied by all the models corre- 
sponding to an envelope with Ti = 4/3,5/3 and 2 respec- 
tively (that is, if the ratio of pressure to energy-density at the 
core-envelope boundary, P b /E b , is set (and not 'arbitrarily' 
assigned) to be about 4.694 x 10~ 2 for all the sequences), the 
corresponding range of the glitch healing parameter turns 
out to be 0.016 < Q < 0.779. This range, however, can ex- 
plain both (the higher as well as lower) values of Q on the 
basis of starquake model of glitch generation for the Crab, as 
well as for the Vela-like pulsars. The sequences correspond- 
ing to an envelope with Ti = 5/3 and 2 yield the maximum 
masses M b ~ 2.06M Q , M c ~ 1.85M© for the Vela (Q ~ 0.2) 
and the minimum masses Mb — 4. IMq, M c ~ 3.98Mq for 
the Crab (Q ~ 0.7) pulsar respectively. The sequence corre- 
sponding to an envelope with Ti = 4/3, however, yields the 
maximum mass M a — 2.67M for the Vela pulsar (Q ~ 0.2) 
but it does not provide the minimum stable mass for the 
Crab pulsar as soon as the constraint of Q > 0.7 is imposed. 

The higher mass values mentioned in the last paragraph 
for the Crab pulsar seem to be unlikely, since none of the ob- 
servational and/or the theoretical study predict such higher 
mass values for the Crab pulsar. Thus, unlike the results 
of steps 2-3 (section 3) which deal with the study of NS 
models without implementing the 'appropriate' fulfillment 
of 'compatibility criterion', the implementation of the 'ap- 
propriate' fulfillment of 'compatibility criterion' also reveals 
that in order to construct a 'realistic' NS sequence composed 
of NS masses comparable with those of the observations, wc 
have to modify the value of matching density at the core- 
envelope boundary. In view of the modern EOSs of dense 
nuclear matter, the upper bound on NS mass compatible 
with causality and dynamical stability can reach up to a 
value as large as 2.2M , since among the variety of mod- 



ern EOSs discussed in the literature only the following EOSs 
yield the maximum mass of NS model in excess of 2M : SLy 
(Douchin & Haensel 2001) EOS, M max = 2.O5M ; BGN1 
(Balberg & Gal 1997) EOS, M max = 2.18M ; and APR (Ak- 
mal et. al. 1998) EOS, M max = 2.21M (see, e.g. Haensel 
et al 2006). However, on the basis of other modern EOSs 
for NS matter, fitted to experimental nucleon-nucleon scat- 
tering data and the properties of light nuclei, Kalogera & 
Baym (1996; and references therein) have also shown that 
the lowest possible upper bound on NS mass, compatible 
with causality and dynamical stability, corresponds to the 
value of 2.2M which can exceed up to a value as large 
as 2.9M . Considering their mean ~ 2.6M as the most 
likely recent value to the upper bound on NS masses and 
substituting this value as an upper bound to our models 
mentioned in step 4 of section 3, we obtain the 'appropriate' 
value of matching density E b = 7.0794 x 10 14 g cm~ 3 at the 
core-envelope boundary. On the basis of this density, our re- 
constructed sequences (corresponding to an envelope with 
Ti = 5/3 and 2 respectively) yield the maximum masses 
M b ~ 1.28M Q , M c ~ 1.14M for the Vela (Q ~ 0.2) and the 
minimum masses Mb ~ 2.53M , M c ~ 2.45M for the Crab 
(Q ~ 0.7) pulsar as shown in Fig. 6. The sequence corre- 
sponding to an envelope with Ti = 4/3, however, yields the 
maximum mass M a ~ 1.66M for the Vela pulsar (Q ~ 0.2) 
but as shown in Fig.6, the minimum mass of the Crab pul- 
sar (M a ~ 2.46M ) belongs to the unstable branch of the 
sequence as soon as the constraint of Q > 0.7 is imposed. 
Obviously, the other conclusions of the present study, drawn 
on the basis of assigning the fiduciary density E b = 2.7 x 10 14 
g cm -3 at the boundary will remain unaltered. 
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Table 1. Mass (M), size (R), compactness ratio (u s tft = M/R) and the central value of the 'local' adiabatic index (ri)o of the 
configuration for various values of the ratio of central pressure to central energy-density (Po/Eo) as obtained by assigning a fiduciary 
value of the density at the surface, E a = 2.7 X 10 14 g cm~ 3 , for the stiffest EOS, P = (E — E s ). It is seen that for each assigned 
value of Po/Eo, the inequality, u st g < (where represents the corresponding value of the compactness ratio for the homogeneous 
density distribution shown in column 2) always remains unsatisfied. The slanted values correspond to the limiting case upto which the 
configuration remains pulsationally 'stable'. 
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Table 2. Mass (M), size (R), compactness ratio (u = M/R), and the 'local' value of adiabatic index at the centre, (ri)o, for different 
values of the ratio of central pressure to central energy-density, (Pq/Eq), for the core-envelope models discussed in the text (steps (2) 
and (3), section 3). Various parameters are obtained by assigning a fiduciary value of Ef, = 2.7 X 10 14 g cm~ 3 for an assigned value of 
the ratio of pressure to energy-density, (Pb/Ei,) = (P^/E^i ~ 1.0645 X 10~ 2 , at the core-envelope boundary. The superscripts a, b, and 
c which appear among various parameters represent the models corresponding to an envelope with Ti = 4/3,5/3, and 2 respectively. 
The slanted values correspond to the limiting case upto which the configurations remain pulsationally 'stable'. It is seen that the M — R 
relation does not provide the necessary and sufficient condition of dynamical stability for any of the sequence, since none of the sequences 
satisfy both of the inequalities, viz. u max < 0.3406 and (ri)o < 2.5946, simultaneously at the maximum value of mass. Although, all 
members of the sequence corresponding to an envelope with Ti = 4/3 fulfill the 'compatibility criterion', u < uj, (where uj, represents 
the compactness ratio of homogeneous density sphere for the corresponding value of Pq/Eq) without fulfilling the necessary and sufficient 
condition of dynamical stability provided by the M — R relation, various members in the stable (as well as unstable) branch of the 
sequences corresponding to an envelope with Ti =5/3 and 2 do not fulfill the 'compatibility criterion' together with the last condition 
provided by the M — R relation as shown in Table 3. Thus, the assigned value of (Pb/Ei,) = (P^/E^i ~ 1.0645 X 10~ 2 corresponds to 
the minimum value for which the 'compatibility criterion', in fact, is not 'appropriately' satisfied by the sequence corresponding to an 
envelope with Ti = 4/3 as discussed in step (3) of section 3. 
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3 


.33138 


17.61968 


0. 


27926 


3. 


,33085 


17.47471 





,28153 


0. 


35897 


3.78575 


3 


.45739 


18.16419 


0. 


28113 


3 


.45653 


17.72636 


0. 


.28801 


3 


,45603 


17.58420 





,29029 


0. 


37791 


3.64616 


3 


.54070 


18.20241 





.28730 


3 


.53990 


17.78229 


0. 


29402 


3 


,53946 


17.65346 





29613 


0. 


43583 


3.29448 


3 


.75240 


18.24864 


0. 


30371 


3 


.75176 


17.88056 


0. 


30991 


3 


,75139 


17.76321 





.31193 


0. 


48898 


3.04509 


3 


.89663 


18.20641 


0. 


31612 


3 


.89610 


17.87641 


0, 


32191 


3 


,89580 


17.76979 





.32381 


0. 


56150 


2.78094 


4 


.02745 


18.03557 





.32982 


4. 


.02703 


17.74399 


0. 


33521 


4. 


,02677 


17.64993 





33697 


0. 


63427 


2.57662 


4 


.09052 


17.74603 





.34045 


4. 


.09018 


17.48404 


0. 


34553 


4. 


08998 


17.40243 





.34713 


0. 


66738 


2.49839 


4 


.09795 


17.57442 


0. 


34440 


4 


.09764 


17.32866 





34926 


4. 


.09745 


17.24674 





.35090 


0. 


68401 


2.46197 


4 


.09650 


17.48041 


0. 


34613 


4. 


.09620 


17.23688 





35100 


4. 


.09604 


17.16103 





.35253 


0. 


70503 


2.41839 


4 


.08963 


17.34875 


0. 


34817 


4. 


.08935 


17.11346 


0. 


35294 


4. 


,08917 


17.03632 





,35452 


0. 


71465 


2.39928 


4 


.08462 


17.28612 


0. 


.34901 


4. 


.08434 


17.05091 


0. 


35380 


4. 


,08419 


16.97883 





35529 


0. 


75963 


2.31644 


4 


.04483 


16.95813 


0. 


35229 


4 


.04458 


16.73528 


0. 


35696 


4. 


,04443 


16.66216 





,35851 


0. 


.81199 


2.23154 


3 


.95964 


16.49323 


0. 


35459 


3 


.95942 


16.28068 


0. 


35920 


3 


,95929 


16.21292 


0, 


,36069 


0. 


83724 


2.19440 


3 


.90077 


16.23405 


0. 


.35490 


3 


,90056 


16.02483 





35951 


3 


90044 


15.95891 





,36099 
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Table 3. Compactness ratio, u(= M/R), and the fractional moment of inertia, Q(= -fcore/Aotal) > f° r different values of the ratio of 
pressure to energy-density, (Pq/Eq), at the centre for the models presented in Table 2. For each value of (Po/Po), the compactness 
ratio of homogeneous density sphere, u^, is shown in column 2. Various parameters are obtained by assigning a fiduciary value of 
Ef, = 2.7 X 10 14 g cm~ 3 for an assigned value of the ratio of pressure to energy-density, (Pj/P;,) = (Pj,/P;,)i ~ 1.0645 X 10~ 2 , at the 
core-envelope boundary. The superscripts a, b, and c which appear among various parameters represent the models corresponding to an 
envelope with Ti = 4/3, 5/3, and 2 respectively. The slanted values correspond to the limiting case upto which the configurations remain 
pulsationally 'stable'. It is seen that for an assigned value of (Po/Po) the inequality, u < u/j, is always satisfied (but not 'appropriately' 
satisfied, cf. the caption of Table 2) by all members of the sequence corresponding to an envelope with Ti = 4/3, whereas it remains 
unsatisfied by various members in the stable (as well as unstable) branch of the sequences corresponding to an envelope with Ti = 5/3 
and 2 respectively. It follows from this table that the choice of P;,/E(,(~ 1.0645 X 10 — 2 ; in the present context, for example) can provide a 
suitable explanation only for the higher values of the glitch healing parameter Q in the range, 0.558 < Q < 0.948 (in the present context, 
for example), on the basis of the starquake mechanism of glitch generation as shown in Fig. 3. 



(P /P ) u h u a u b u c Q a Q b Q c 



0. 


10603 





14794 





12862 





.14020 


0. 


14389 





.55820 





.65183 


0. 


.68421 





11977 





.16070 


0. 


.14250 





.15367 


0. 


15725 





.60171 





.68740 





.71684 





12544 





.16567 


0. 


14789 


0. 


15897 





16241 





61720 





.70036 





.72805 


0. 


13235 





17152 


0. 


15435 


0. 


16514 





16860 





.63523 





.71436 





.74134 


0. 


15071 


0. 


.18603 


0. 


.17027 





.18057 


0. 


18390 





.67538 





.74645 


0. 


.77083 


0. 


18604 


0. 


21029 


0. 


.19699 


0. 


.20654 


0. 


20959 





.73256 





.79226 


0. 


.81213 


0. 


.22012 





23000 


0. 


21889 


0. 


.22770 





23056 





77223 





.82325 





.84041 





25219 





.24591 


0. 


23655 


0. 


24486 





24756 





80003 





.84534 





.86049 


0. 


28047 





25823 





25025 


0. 


.25817 


0. 


26065 


0. 


.81961 





.86085 


0. 


.87416 





.31322 





.27081 


0. 


26427 


0. 


.27172 


0. 


.27416 





.83806 





.87524 


0. 


.88769 





33333 





.27778 


0. 


27209 


0. 


.27926 





.28153 





84795 





.88286 





.89418 





35897 





28593 


0. 


.28113 


0. 


.28801 





29029 





.85858 





.89113 





.90219 


0. 


37791 





29149 


0. 


.28730 





.29402 


0. 


29613 





.86561 





.89685 


0. 


.90686 


0. 


.43583 





30640 


0. 


30371 


0. 


.30991 


0. 


31193 


0. 


.88305 





.91052 


0. 


.91962 





48898 





31785 


0. 


.31612 





.32191 


0. 


.32381 





.89536 





.92014 


0. 


.92842 





56150 





33083 


0. 


32982 


0. 


.33521 





.33697 


0. 


.90837 





.93055 





.93793 


0. 


.63427 





.34152 


0. 


34045 


0. 


.34553 


0. 


.34713 





.91811 





.93842 





.94493 


0. 


66738 





34577 





.34440 





.34926 





.35090 





.92180 





.94106 





.94765 


0. 


68401 





34778 


0. 


.34613 


0. 


.35100 


0. 


35253 





.92331 





.94251 


0. 


.94864 





70503 





35021 


0. 


34817 


0. 


.35294 





35452 


0. 


.92529 





.94399 





.95028 





.71465 





.35128 


0. 


34901 


0. 


35380 





35529 


0. 


.92604 





.94480 





.95070 


0. 


.75963 





35600 


0. 


.35229 





.35696 


0. 


35851 





.92925 





.94739 


0. 


.95349 


0. 


81199 





36095 


0. 


35459 





.35920 


0. 


36069 





.93215 





.94997 


0. 


.95578 





83724 





.36314 


0. 


35490 


0. 


.35951 





36099 





93292 





.95075 





.95650 
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Table 4. Mass (M), size (R), compactness ratio (u = M/R), and the 'local' value of adiabatic index at the centre, (ri)o, for different 
values of the ratio of central pressure to central energy-density, (Pq/Eo), for the core-envelope models discussed in the text (step (4), 
section 3). Various parameters are obtained by assigning a fiduciary value of E b = 2.7 X 10 14 g cm -3 and for an 'appropriate' value of 
the ratio of pressure to energy-density, (P b /E b ) = {Pb/^b)2 — 4.694 X 10~ 2 , at the core-envelope boundary. The superscripts a, b, and 
c which appear among various parameters represent the models corresponding to an envelope with Ti = 4/3,5/3, and 2 respectively. 
The slanted values correspond to the limiting case upto which the configurations remain pulsationally stable. It is seen that the M — R 
relation does provide the necessary and sufficient condition of dynamical stability for all of the sequences, since all of the sequences 
satisfy both of the inequalities, viz. u max < 0.3406 and (ri)o < 2.5946, simultaneously at the maximum value of mass. Because all 
members of the sequences corresponding to an envelope with Ti = 4/3,5/3 and 2 fulfill the 'compatibility criterion', u < uj, (where 
«h represents the compactness ratio of homogeneous density sphere for the corresponding value of Pq/Eq), together with fulfilling the 
necessary and sufficient condition of dynamical stability provided by the M — R relation, it follows, therefore, that the assigned value of 
(P b /E b ) = {P b /E b ) 2 ~ 4.694 X 10~ 2 corresponds to the minimum value for which the 'compatibility criterion' is 'appropriately' satisfied 
by all the sequences corresponding to NS models with an envelope Ti = 4/3, 5/3 and 2 respectively as discussed in step (4) of section 3. 



(Po/E ) 


(Ti)o 


(M«/M ) 


R a (km) 


u a 


(M b /M Q ) 


R b (km) 




u» 


(M c /M ) 


R c {km) 




u c 


0. 


10052 


10.9487 


1. 


.79122 


43 


.32919 


0.06106 


1. 


55535 


20.66582 





11116 


1. 


.47601 


17. 


50211 





.12456 


0. 


11236 


9.90023 


1. 


.87195 


39 


.26268 


0.07042 


1. 


67909 


20.40422 


0. 


.12154 


1 


61142 


17. 


58182 


0. 


13537 


0. 


12029 


9.31353 


1. 


93294 


37 


11014 


0.07693 


1. 


76262 


20.27094 





.12843 


1 


.70131 


17. 


65121 


0. 


.14236 


0. 


.12955 


8.71910 


2. 


00865 


35 


.04138 


0.084665 


1. 


85989 


20.15492 





.13630 


1 


.80486 


17. 


.73130 





.15034 


0. 


.14549 


7.87332 


2. 


14512 


32 


.33658 


0.09798 


2. 


02462 


20.01554 





.14940 


1 


.97821 


17. 


.88177 


0. 


.16340 


0. 


15354 


7.51305 


2. 


21518 


31 


.27768 


0.10461 


2. 


10582 


19.97345 





.15572 


2 


06296 


17. 


95627 


0. 


.16969 


0. 


20061 


5.98487 


2. 


61355 


27 


42286 


0.14077 


2. 


54549 


19.87280 


0. 


.18919 


2 


.51677 


18. 


35244 


0. 


.20255 


0. 


23990 


5.16877 


2. 


.91143 


25 


.75275 


0.16698 


2. 


86158 


19.87300 





.21268 


2 


.83972 


18. 


61125 


0. 


.22536 


0. 


27901 


4.58413 


3. 


.17036 


24 


.69328 


0.18963 


3. 


13199 


19.87624 





.23274 


3 


.11471 


18. 


80100 


0. 


.24469 


0. 


.34816 


3.87220 


3. 


54016 


23 


.49470 


0.22255 


3. 


51403 


19.84340 





.26156 


3 


.50187 


18. 


97877 


0. 


.27253 





44840 


3.23013 


3. 


90788 


22 


.35809 


0.25816 


3. 


89112 


19.64414 





29256 


3 


.88312 


18. 


97775 


0. 


30222 


0. 


48306 


3.07012 


3. 


99636 


22 


.02864 


0.26795 


3. 


98169 


19.54181 





30094 


3. 


.97464 


18. 


92270 





.31024 


0. 


51773 


2.93149 


4. 


06762 


21 


.71508 


0.27667 


4. 


05467 


19.41360 





30848 


4. 


.04842 


18. 


.83572 


0. 


.31746 


0. 


.54455 


2.83637 


4. 


11157 


21 


.47901 


0.28273 


4. 


09975 


19.30104 





.31373 


4. 


.09402 


18. 


74922 





.32251 


0. 


58709 


2.70330 


4. 


16243 


21 


.10649 


0.29128 


4. 


15213 


19.09782 





32112 


4. 


14711 


18. 


58635 


0. 


32956 


0. 


63150 


2.58354 


4. 


19170 


20 


.70981 


0.29895 


4. 


18272 


18.84735 





.32778 


4. 


.17831 


18. 


37295 


0. 


.33589 


0. 


66866 


2.49552 


4. 


.19788 


20 


.36982 


0.30438 


4. 


18981 


18.61628 





33242 


4 


,18584 


18. 


.16228 





.34040 





.68296 


2.46421 


4. 


19565 


20 


.23663 


0.30623 


4. 


18789 


18.51726 





.33404 


4. 


18408 


18. 


07316 


0. 


.34194 


0. 


69900 


2.43062 


4. 


19032 


20 


.08121 


0.30820 


4. 


18289 


18.40185 





.33573 


4. 


.17923 


17. 


96912 


0. 


.34352 


0. 


.72302 


2.38309 


4. 


17593 


19 


.84795 


0.31076 


4. 


16896 


18.22159 





.33793 


4. 


16552 


17. 


79888 


0. 


.34567 


0. 


74010 


2.35117 


4. 


16090 


19 


.67543 


0.312351 


4. 


15423 


18.08144 





.33934 


4. 


.15094 


17. 


66827 





.34700 





77503 


2.29027 


4. 


11752 


19 


.30611 


0.31501 


4. 


11141 


17.77494 





34164 


4. 


.10839 


17. 


38013 





.34914 


0. 


80605 


2.24061 


4. 


06269 


18 


.95854 


0.31651 


4. 


05700 


17.47310 





.34294 


4 


.05419 


17. 


08740 


0. 


35044 


0. 


.83972 


2.19088 


3. 


98260 


18 


.54954 


0.31711 


3. 


97730 


17.10352 





34347 


3 


.97469 


16. 


72970 


0. 


35091 
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Table 5. Compactness ratio, u(= M/R), and the fractional moment of inertia, Q(= -fcore/Aotal) > f° r different values of the ratio of 
pressure to energy-density, (Pq/Eq), at the centre for the models presented in Table 4. For each value of (P0/-E0), the compactness 
ratio of homogeneous density sphere, u^, is shown in column 2. Various parameters are obtained by assigning a fiduciary value of 
Ef, = 2.7 X 10 14 g cm -3 and for an 'appropriate' value of the ratio of pressure to energy-density, (Pb/Ei,) = {P\,/E} > )2 — 4.694 X 10 -2 , 
at the core-envelope boundary. The superscripts a, b, and c which appear among various parameters represent the models corresponding 
to an envelope with Y\ = 4/3,5/3, and 2 respectively. The slanted values correspond to the limiting case upto which the configurations 
remain pulsationally stable. It is seen that for an assigned value of (Pq/Eq), the inequality, u < u^, is 'appropriately' satisfied by 
all members of the sequences corresponding to an envelope with Ti = 4/3,5/3 and 2 respectively. It follows from this table that the 
'appropriate' choice of Pj,/.E(,(~ 4.694 X 10 ; in the present context, for example) can provide a suitable explanation for both (the 
higher as well as lower) values of the glitch healing parameter Q in the range, 0.016 < Q < 0.779 (in the present context, for example), 
on the basis of the starquakc mechanism of glitch generation as shown in Fig. 5. 



(Po/Eo) u h u a u b u c Q a Q b Q c 



0. 


10052 





14253 


0.06106 


0. 


11116 





.12456 


0. 


01548 


0. 


06873 





.09775 


0. 


11236 





.15394 


0.07042 


0. 


.12154 





.13537 





02739 


0. 


09920 


0. 


.13475 


0. 


12029 





16116 


0.07693 


0. 


.12843 


0. 


14236 


0. 


.03739 


0. 


12071 


0. 


15969 


0. 


12955 





16918 


0.08466 





13630 


0. 


15034 





.05089 


0. 


14626 


0. 


18860 


0. 


.14549 





.18205 


0.09798 


0. 


14940 





16340 





07777 


0. 


19020 





.23645 


0. 


15354 





18814 


0.10461 


0. 


15572 





16969 


0. 


09261 


0. 


21182 





.25950 





20061 





.21911 


0.14077 





18919 


0. 


20255 





18555 


0. 


32652 


0. 


.37684 





23990 





.24008 


0.16698 


0. 


.21268 





22536 





26028 


0. 


40470 


0. 


.45366 





.27901 





.25763 


0.18963 


0. 


.23274 





24469 





32671 


0. 


46877 





.51517 





34816 


0. 


.28259 


0.22255 


0. 


.26156 





27253 





.42326 


0. 


55590 





.59803 


0. 


44840 





.30929 


0.25816 


0. 


.29256 


0. 


30222 





.52531 


0. 


64342 


0. 


.67967 





48306 





.31666 


0.26795 





.30094 


0. 


31024 





.55283 


0. 


66607 


0. 


.70070 


0. 


51773 





32332 


0.27667 


0. 


30848 





.31746 





57715 


0. 


68625 





.71944 


0. 


.54455 





32803 


0.28273 


0. 


.31373 


0. 


32251 





59401 


0. 


70020 





.73250 


0. 


58709 





.33482 


0.29128 





.32112 


0. 


32956 





.61777 


0. 


71982 


0. 


.75062 


0. 


63150 


0. 


.34115 


0.29895 


0. 


.32778 


0. 


33589 


0. 


.63925 


0. 


73770 


0. 


.76712 





66866 





34593 


0.30438 





.33242 


0. 


34040 





.65474 





.75034 





77916 


0. 


68296 





.34766 


0.30623 





.33404 


0. 


.34194 





.66008 


0. 


75491 


0. 


.78336 


0. 


69900 


0. 


34952 


0.30820 


0. 


.33573 


0. 


34352 


0. 


66592 


0. 


75976 


0. 


.78778 





72302 





35220 


0.31076 


0. 


.33793 


0. 


34567 





.67364 


0. 


76625 


0. 


.79403 





.74010 





.35401 


0.31235 


0. 


33934 





34700 





67868 


0. 


77068 


0. 


.79815 


0. 


77503 





.35751 


0.31501 


0. 


.34164 





34914 





68772 


0. 


77855 





.80543 





80605 





.36041 


0.31651 





.34294 


0. 


35044 





.69398 


0. 


78427 


0. 


.81111 


0. 


.83972 





.36336 


0.31711 


0. 


.34347 


0. 


35091 


0. 


69887 


0. 


78917 


0. 


.81587 



